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The thermoelectric properties of quantum dot arrays (QDAs) embedded in nanowires connected to 
electrodes are studied theoretically in the Coulomb blockade regime. A Hurbbard-Anderson model 
is used to simulate the electronic contribution to thermoelectric properties of a QDA junction 
system. The electrical conductance, Seebeck coefficient, and electron thermal conductance are 
calculated by both the Keldysh Green function method and the mean-field approach. The phonon 
thermal conductivities are calculated by using the equation of phonon radiative transfer method. In 
the Coulomb blockade regime the electron thermal conductance is much smaller than the phonon 
thermal conductance. Therefore, the optimal figure of merit (ZT) can be enhanced by increasing 
thermal power and decreasing phonon thermal conductance simultaneously. We found that it is 
possible to obtain ZT value of InGaAs/GaAs QDAs embedded in nanowires larger than one at 
room temperature. 


I. INTRODUCTION 

Recently, many efforts have been devoted to seeking 
efficient thermoelectric (TE) materials with the figure 
of merit {ZT) larger than 3, because there are poten¬ 
tial applications of solid state thermal devices such as 
coolers and power generators, which can replace con¬ 
ventional compressor-based refrigerators and fossil fuel 
generators to reduce CO 2 emission [1-7]. Nevertheless, 
the optimization of TE materials is extremely difficult, 
since ZT = S^GeT/n depends on the electrical conduc¬ 
tance (Ge), Seebeck coefficient (S),and thermal conduc¬ 
tance {k). T is the equilibrium temperature. These 
physical quantities are usually related to one another. 
Mechanisms leading to the enhancement of power factor 
{PF = S^Ge) would also enhance the thermal conduc¬ 
tance. Consequently, it is difficult to obtain ZT above 
one in conventional bulk materials. [1] 

Impressive ZT values for quantum dot array (QDA) 
embedded in nanowires have been experimentally demon¬ 
strated. [8] The power factor and thermal conductance 
become independent thermoelectric variables under the 
condition Ke/uLa “C 1, where Ke and KLa denote, re¬ 
spectively, the electron thermal conductance and lat¬ 
tice thermal conductance. [8] In the Coulomb blockade 
regime, electron transport process is seriously suppressed 
by the electron Coulomb interactions, therefore Ke as well 
as Ge are reduced significantly. [9] Under the condition 
i^e/uLa ^ 1, one can increase the power factor and de¬ 
crease the phonon thermal conductance simultaneously 
to optimize 2'T.[9] 

Thermoelectric properties of quantum dots (QDs) em¬ 
bedded in a matrix connected to metallic electrodes were 
studied by several groups in the absence of phonon ther¬ 
mal conductivity. [10-16] For the applications of solid 
state coolers and power generators at room temperature, 
one needs to consider a large number of serially coupled 


QDs, otherwise it is not easy to maintain a large tempera¬ 
ture difference across the QD junction, which was pointed 
out to be crucial in the implementation of high-efRciency 
thermoelectric devices. [1,2] In addition, the phonon ther¬ 
mal conductivity plays a significant role in the optimiza¬ 
tion of ZT at high temperatures. In this paper, we carry 
out theoretical analysis of ZT of QD arrays embedded in 
nanowires, including the phonon conductance, which is 
calculated by using the phonon radiative transfer method 
as introduced in Ref. [17]. Although the method does not 
take into account the microscopic mechanisms associated 
with quantum confinement of acoustic phonons, it gives 
reasonable agreement with the lattice-dynamics model 
and experiments for nanowires by merely considering the 
boundary scattering effect of phonons. 


II. FORMALISM 

A QDA embedded in a nanowire connected to the 
metallic electrodes can be described by the Hurbbard- 
Anderson model. Here we consider nanoscale semicon¬ 
ductor QDs, in which the energy level separations are 
much larger than their on-site Coulomb interactions and 
thermal energies. Thus, only one energy level for each 
quantum dot needs to be considered. The Hamiltonian 
of the system is given by 77 = iJo + Hqd'- 

Ho = +'^ekbl^^bk,a (1) 

k,(T k,(T 

fe,(7 fc,(7 

where the first two terms describe the free electron gas of 
left and right electrodes. a\ ^ {b\ „) creates an electron of 
momentum k and spin a with energy Ck in the left (right) 
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electrode. 14/ [I = L, R) describes the coupling between 
the electrodes and the left (right) QD. d] ^ [di,a) creates 
(destroys) an electron in the .^-th dot. 


and 

^22 = ^1 deT{e){e-Ep){^)E,. 


( 10 ) 


Hqd = ^ ^ UiUt^ane^a ( 2 ) 

i,a a 

l,j,CT,CT' 

where Ee is the spin-independent QD energy level, and 
ni^cr = dg g.df,cr- Ui and Uej describe the intradot and in¬ 
terdot Coulomb interactions, respectively, tej describes 
the electron interdot hopping. Note that the interdot 
Coulomb interactions as well as intradot Coulomb inter¬ 
actions play a significant role on the charge transport for 
semiconductor QD array. 

Using the Keldysh-Green’s function technique [18], the 
charge and heat currents of electrons leaving electrodes 
are expressed as 


J 

Ql{r) 


^ I der(e)[/i(e)-/«(e)], (3) 

dL{R))[fL{e) - fR{e)],{4:) 


where T(e) is the transmission coefficient. fL(R){d) = 
+ 1] denotes the Fermi distribution 
function for the left (right) electrode, and hr denote 
the chemical potentials of the left and right leads, respec¬ 
tively, with their average denoted by Er = {hl + Hr)I‘^- 
{hl — Hr) = eAV is the voltage across the QDA junc¬ 
tion. Tl(_r) denotes the equilibrium temperature of the 
left (right) electrode, e and h denote the electron charge 
and Planck’s constant, respectively. Ql{r) denotes the 
heat current leaving from the left (right) electrode. 

In the linear response regime, Eqs. (3) and (4) can be 
rewritten as 



. at 

(5) 


+ ^ 12 ^ 
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where there are two sources of driving force to yield the 
charge and heat currents. AT = Tl — Tr is the temper¬ 
ature difference across the junction. The thermoelectric 
coefficients in Eqs. (5) and (6) (Tii, T 12 , -C 21 , and £ 22 ) 
are evaluated by 


£11 — 


£12 = 


2e‘^T 

h 

2eT^ 


I dene){^)E,, ( 8 ) 


£21 = ^ / *r(e)(e - Ee){^)t, (9) 


Here T(e) and /(e) = l/[e('^“^^)ABT _|_ evaluated 

under the equilibrium condition. The detailed expression 
of T(e) can be found in Ref. [19]. 

If the system is in an open circuit, the electrochemical 
potential will be established in response to a temperature 
gradient; this electrochemical potential is known as the 
Seebeck voltage (Seebeck effect). The Seebeck coefficient 
(amount of voltage generated per unit temperature gra¬ 
dient) is defined as 5 = AV/AT = —£i 2 /(T£ii). To 
judge whether the system is able to generate or extract 
heat efficiently, we need to consider the figure of merit 

|i| 

= (zr)„ 

Ke+KLa ^+KLa/Ke 

Here, Ge = EwjT is the electrical conductance and 
Ke = ((£ 22 /'£^) — S'^GeT) is the electron thermal con¬ 
ductance. {ZT)q represents the ZT value in the absence 
of phonon thermal conductance, klq. In the Hamilto¬ 
nian H, the term [Ele^ph) describing the interactions be¬ 
tween electrons and phonons is ignored. For simplicity, 
we adopt KLa = RwireFa to describe the phonon ther¬ 
mal conductance of QDA embedded a nanowire. The 
dimensionless scattering factor i4 is used to include the 
phonon scattering effect arising from surface boundary of 
QDs.[l] It is possible to reduce phonon thermal conduc¬ 
tance by one order of magnitude when QDs in a quan¬ 
tum wire behave like phonon scatterers.[20] Therefore, 
the maximum Eg is assumed to be 0.1 in this study. For 
example, Ref. [20] pointed out that the phonon ther¬ 
mal conductance reduction of QDA nanowires in a high 
temperature regime arises from the filtering of high fre¬ 
quency phonons. [2] However, the method considered in 
Ref. [20] requires a heavily numerical calculation of 
in the optimization of ZT. 

Based on Holland’s model [21,22] in which the longi¬ 
tudinal and transverse acoustical phonon branches are 
treated separately, we calculate K^jire of a rectangular 
nanowire with cross sectional area A and length Lx- The 
thermal conductance Kwire is related to the phonon ther¬ 
mal conductivity by Hwire = RphA/Lx^ where Kph = 
Ktq + Ktu + Kr denotes the thermal conductivity with 


Kt„ = -T 
° 3 


Grox'^e^ie^ -l)-^dx 


+Tj^ + PtxT^ ’ 


( 12 ) 
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,(13) 


and 

= 1^3 CLx'^e-je- - l)-^dx 


(14) 
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The phonon thermal conductivities Ktqi Ktu, and 
result from the low-frequency transverse branch, high- 
frequency transverse branch, and longitudinal acous¬ 
tical phonon branch, respectively. In Eqs. (12)-(14), 
we have Cj=To,Tu,L = {ks/K)^, where Vj 
is the phonon group velocity for the j*^-branch. 
is the phonon boundary scattering rate, and = 

ax'^T'^ is the phonon-impurity scattering rate. 
j3TuX^T'^/ sinh(a;), and Plx^T^ arise from three-phonon 
scattering rates, a and f3j are empirical parameters de¬ 
termined by fitting experimental data. 9i{TujJi/kB){i = 
1,2, and 3) denotes the Debye temperature. The fre¬ 
quencies of LOi, for Tq, Tu and L modes can be found in 
Ref. [17]. The a, /3j, Vj and 9i parameters for Si and 
GaAs are adopted from Refs. [17] and [22], respectively. 
Eqs. (12)-(14) can well describe the phonon thermal con¬ 
ductivity of semiconductor materials with Tb = Vb/L. 
Vb{= [[2v ^^is the average phonon group ve¬ 
locity, where vt and vl correspond to the group velocities 
of the transverse and longitudinal branches. L denotes 
the sample size. Chen and Tien [21] extended the Holland 
model to illustrate the phonon thermal conductivity of 
quantum wells by considering the geometry effect on r^. 
Based on the formalism of ref. [21], the phonon-boundary 
scattering rate of quantum wire is derived and de¬ 

termined by 


= (15) 


where Tt = -I- (Hi -|- B 2 )T^uj^ + VblLc)~^ is the 

total internal relaxation time of the bulk which includes 
normal process and umklapp process. A = a(fi/kB)^, 
and Bj = /3j{h/kB)^- Lc denotes the sample length. 
The factor of Q = -I- Q ^ + Q illustrates 

the boundary effects on the phonon thermal conductivity 
of nanowires. The expression of Q is given by 
(i) 0 < 0 < 7r/2 and 0 < (j) < n 


g++/u 


dfl 


(e — l)(e — 1) 
— 1 ’ 


(16) 


(ii) 0 < 9 < Ti 12 and tt < (j) < 2 tt 


g+-/n 


dn 


(e 


— 


l)(e-«^ -1) 
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(17) 


(hi) TT j2 < 9 < TT and 0 < (f> < tt 


g-+/n 


dn 


(e 


— 


l)(e-«^ - 1) 


3 Ce — 1 


(18) 


(iv) TT j2 < 9 < TT and tt < (ft <2 tt 


g—/A 


— 1 


(19) 


We have the notations H = ^ and dfl = 

dydzd9d(l) sin^(0) cos'^{(j)) in Eqs. (16)-(19), where Ly(^z) 


is the lateral size of rectangular wire in the y(z) direc¬ 
tion. The other notations are = —y/(A sin(0) sin(()))), 

= -2:/(Acos(6>)), ^y = {Ly - ?/)/(Asin(6»)sin(()))) and 
^z = {Lz — z)/{Kcos{9)). The average phonon mean free 
path is assumed to be A = VbTt- 

The phonon thermal conductivities {Kph) of rectangu¬ 
lar Si and GaAs nanowires calculated by using Eqs. (12)- 
(14) are plotted for two different topological structures 
in Fig. 1. The calculations of silicon nanowire can be 
used to compare with experimental results to examine 
the validity of the scheme adopted in Ref. [21]. The heat 
problem of high efficiency solar cells made of III-V com¬ 
pounds (such as GaAs) can degrade the system perfor¬ 
mance. [1] Thus, it is important to design a GaAs cooler, 
which can be integrated with the solar cell to solve such 
a problem. Figure 1(a) shows Kph of Si nanowire with 
Ly = Lz as a function of temperature for three cross- 
section areas: solid line (A = (30 nmf), dashed-dotted 
line (A = (10 nmY) and dashed line (A = (5 nm)'^). 
It is expected that Kph decreases with decreasing cross- 
section area due to enhanced boundary scattering. For 
a nanowire with small cross section, Kph increases slower 
with respect to kBT for T > 200 K. These results are 
consistent with other theoretical calculations [23] and ex¬ 
perimental observations.[24,25] For nanowires with large 
cross section, however, Kph is slightly underestimated 
compared to other theoretical works.[23] For the smallest 
area considered (A = (5 nm)^), the Kph value obtained 
is very close to the result calculated from the lattice- 
dynamics model in Ref. 20, which takes into account the 
quantum confinement of acoustic phonons. Although a 
nanowire with even smaller cross section leads to much 
smaller phonon thermal conductivity, it is a challenge to 
implement such tiny nanowires.[1,2] 

Figure 1(b) shows Kph of GaAs nanowire as a func¬ 
tion of temperature for different L^ sizes with Ly fixed at 
5 nm. We note that Kph of GaAs is smaller than that of Si 
in a wide temperature range for A = (5 nmY because the 
average group velocity in GaAs is smaller. Kph increases 
with increasing L^. However, Kph for Lz = 500 nm be¬ 
comes almost the same as Kph of GaAs thin film with 
Lz = Lc = 0 . 729 cto and Ly = b nm (see the solid and 
black line). This implies that the size effect of Lz on Kph 
can be ignored when Lz is larger than 500 nm, compara¬ 
ble to the average phonon mean free path of GaAs. The 
maximum Kph value of GaAs thin film at T = lOO AT is 
around 4 WjmK, which is still smaller than that of Si 
nanowire with A = (10 nm)^. This implies that GaAs 
may have a better ZT value than silicon at room tem¬ 
perature. Our results for GaAs are consistent with the 
calculation of Ref. [26]. Although the TE properties of 
germanium/silicon QDs system are interesting for their 
low cost fabrication process [27], the silicon semiconduc¬ 
tor TE devices can not be directly integrated with HI-V 
compound solar cell systems. [1] Therefore, we focus on 
the ZT optimization of InGaAs/GaAs QDs junction sys¬ 
tem in the next section. 
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III. RESULTS AND DISCUSSION 

Because it is important to reduce the phonon ther¬ 
mal conductivity in the optimization of ZT, we need 
to consider a QDA nanowire with long length. There¬ 
fore, before investigating the TE properties of QDAs with 
large number of coupled dots, we first consider the sim¬ 
plest structure of QDA, which consists of serially cou¬ 
pled triple QDs (SCTQD). Although we have theoret¬ 
ically studied the transport and TE properties of triple 
QD molecules with full solution in the Coulomb blockade 
regime,[28,29] it is still a challenging work to obtain the 
full solution to the charge transport through QDAs with 
more than three dots due to the fact that the calculation 
scales up exponentially. To solve such a difficult problem 
many efforts which adopt some approximation schemes 
were attempted.[9,30-32] The Hartree-Fock approxima¬ 
tion (HFA) is often used to investigate the TE coeffi¬ 
cients of realistic molecules. [30-32] As a consequence, it 
is desirable to examine the validity of such a mean-field 
approach. At high temperatures, in the absence of in¬ 
terdot Coulomb interactions the results based on the ap¬ 
proach of Ref. [9] can achieve good agreement with those 
calculated by the full solution as reported in Ref. [29]. 
Therefore, we can examine the validity of HFA method 
by using the approach of Ref. [9] instead of using the full 
solution of Ref. [29], which requires much heavier numer¬ 
ical calculation. In addition, the approximation of Ref. 
[9] can provide a closed form transmission coefficient of 
Eqs. (7)-(10), which is useful for analyzing the charge 
transport properties. 

Because the optimization of power factor (PE) fa¬ 
vors identical QDs, size-independent electron hopping 
strength, and symmetrical tunneling rates, [9] in Fig. 2 we 
plot Ge, S', Ke and PF of SCTQD as functions of gate 
voltage {eVg) at two different temperatures by consider¬ 
ing = tc = STq, Ui = Uo = 30ro and symmetrical 
tunneling rates (T^ = T^ = ITo). All energy scales are 
in units of Tq, which is taken to be 1 meV. Here, the 
interdot Coulomb interactions are ignored for simplicity 
when we examine the validity of HFA method. The solid 
lines and dot-dashed lines are calculated according to the 
procedure of Ref. [9] and the HFA method of Ref. [30-32], 
respectively. The approach of Ref. [9], which neglects 
the interdot Coulomb interactions, is called the on-site 
Coulomb approximation (OCA). Both OCA and HFA are 
efficient tools for studying the charge transport through 
molecular junctions. 

For all four physical quantities (Ge, S, Ke and PF = 
S^Ge), the spectra consist two main groups of features, 
which are mirror image of each other with respect to 
the mid point {eVg = dSTo). This is a consequence of 
the electron-hole symmetry possessed by the SCTQD 
structures in the absence of interdot Coulomb interac¬ 
tions [29]. At low temperature (fc^T = O.lPo), there is 
a gap between the two groups, which is related to on¬ 
site Coulomb interaction, Uq. The first three peaks of 
Ge (shown by the solid line) correspond to three reso¬ 


nant channels at Eq — V2tc, Eq, and Eq + \/2tc. For the 
higher-energy group, there are also three peaks arising 
from resonant channels at Eq + Uq — v^tc, Eq + Uq and 
Ec + Uq + V^tc- We find that the behaviors of Ge, S 
and Ke calculated by HFA (dot-dashed lines) are quite 
different from those calculated by OCA (solid lines). 
For example, the maximum Ge of dot-dashed lines is 
larger than e^/h (quantum conductance), but not for 
OCA. This is due to the inability of HFA to include 
electron correlation. The energy positions of resonant 
channels within HFA depend on the average occupation 
number in each QD.[30-32], which can be a fractional 
number. In contrast, within OCA the energy positions 
of resonant channels are determined by integer charges, 
while electron correlation functions and fractional occu¬ 
pation numbers only appear in the probability weights of 
quantum paths. [9] In addition, we observe an unphysi¬ 
cal enhancement of S in the middle of the Coulomb gap 
{eVg = dSTo) for HFA results (dot-dashed line). This 
illustrates the drawback of HFA calculation for TE co¬ 
efficients of QDs at low temperatures when the charging 
process plays a significant role. At high temperatures 
(for example, = lOTo), the detailed resonant struc¬ 
tures are washed out, and the unphysical enhancement 
of S in the middle of the Coulomb gap also disappears. 
These two approaches give qualitatively similar results at 
high temperatures. However, the HFA approach tends 
to overestimate Ge, Ke, and PF. If we include the in¬ 
terdot Coulomb interactions, the interdot correlation ef¬ 
fects become crucial [29]. The symmetrical behavior of 
TE coefficients will be lost. In general, we will obtain 
smaller PF for QD energy levels below Ep than that for 
QD energy levels above Ff’.[9,29] This indicates that to 
achieve maximum PF, it is preferable to have the or¬ 
bital depletion situation with the total occupation num¬ 
ber N = + Nc,a + < 1, where Ni^a- de¬ 

notes the single particle occupation number in dot ^.[29] 

For TE devices operated at room temperature, it is 
important to optimize their ZT values at high temper¬ 
atures. To further clarify the differences between these 
two approaches. Fig. 3 shows Ge, S, Ke and PF as func¬ 
tions of the detuning energy {A = Eq — Ep > 0) at 
various temperatures. Note that in Fig. 3, we consider 
the case with QD energy levels above Ep, which satis¬ 
fies the condition for orbital depletion. In one-particle 
transport process (QD molecule with an empty state) 
only the intradot Coulomb interaction effects are impor¬ 
tant. [29] This is a typical feature for carrier transport 
in the Coulomb blockade regime. [9] The other physical 
parameters are the same as those of Fig. 2. From the 
results of S and PE, these two approaches agree well 
with each other when QD energy levels are far above 
Ep. This is expected, because the electron Coulomb in¬ 
teractions become unimportant for large A values. The 
maximum PF value decreases with increasing tempera¬ 
ture. For kpT = ITo, ksT = STo and kpT = 13ro, 
the maximum PF occurs at near A = STo, A = ISTo, 
and A = 35ro, respectively. When kpT = 26ro, the 
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maximum PF occurs at near A = SOFo- 

From results of Fig. 3, we see that HFA is a reasonably 
good approximation for analyzing the TE coefficients of 
QDAs with QD energy levels above Ep. Thus, we can 
use HFA to investigate TE coefficients of QDA nanowire 
with large numbmer of dots. Here, we use HEA (in¬ 
stead of OCA) of Ref. [30] to study QDA nanowire with 
N=10 because of its simplicity. Fig. 4 shows the Ge, S 
and {ZT)q as functions of temperature for various de¬ 
tuning energies (A). The other physical parameters are 
the same as those of Fig. 2. The behavior of Ge with 
respective to temperature can be understood by consid¬ 
ering Ge ~ l/iksT cosh‘s{A/2kBT) for Ue = 0 (see the 
curve with triangle marks). S is proportional to T when 
temperature approaches zero (not shown here), whereas 
S ~ A/kpT at high temperatures and A/F 1 (see the 
curves with triangle marks). In the absence of Kpa, a 
very large {ZT)q appears in the low temperature regime. 
The maximum (ZT)o is beyond three, which is a critical 
value for the realization of solid state coolers and power 
generators. {ZT)o larger than 3 is a direct consequence of 
the fact Ge ^ Ke {kc not shown here) and large enhance¬ 
ment of 5^ in the Coulomb blockade regime. Using small 
nanowires in the quantum thermal conductance regime 
and a vacuum space to blockade surrounding phonon heat 
current was proposed to yield a vanishingly small phonon 
thermal conductivity. [33] However, such a design requires 
advanced nanotechnology. The results of Fig. 4(c) also 
imply that the on-site electron Coulomb interactions can 
be ignored only when U^/ksT » 1. {ZT)o is highly 
overestimated at high temperatures in the absence of Ue. 

Fig. 5(a) shows ZT in the presence of Kpa (klo = 
FsKph{A/Lx)) for the case of a GaAs nanowire with Lz = 
Ly = 5 nm and Lx = 250 nm. Eg is assumed to be 0.1 
in the calculation. The curves of Fig. 5(a) have one-to- 
one correspondence to those of Fig. 4(c). The maximum 
ZT value is significantly reduced in the presence of Kpa- 
This reduction of ZT results from the fact klo ^ Ke 
[see Eq. (11)], which implies that the optimization of ZT 
can be achieved by reducing Kph and increasing power 
factor PF = S^Ge simultaneously. The maximum ZT 
depends on the detuning energy, A. This behavior can be 
understood from the results of Fig. 3(d). The ZT value is 
almost independent of temperature for A = 40Fo, when 
kpT > 20ro. This is because the ratio of Ge/iTupa) is 
nearly constant. In the high temperature regime we have 
Kph ~ 1/U and Ge decays very slowly with increasing 
temperature. 

Fig. 5(b) shows ZT as a function of temperature for 
various electron hopping strengths at A = dOFo. ZT 
is enhanced with increasing tc- Under the condition 
Kpa S> Kej the behavior of ZT with respect to G can 
be analyzed by the power factor PF = S^Ge- When the 
QD energy is far above Ep (A = dOFo), the behavior of 
S = A/kpT is independent of tc in the Coulomb blockade 
regime once kpT > lOFo. Therefore, the enhancement 
of ZT is completely determined by the increase of Ge 
with respect to tc- When tc is smaller than the tunneling 


rate (Fi, = F/j = F), Ge increases quickly with increas¬ 
ing tc, and it becomes insensitive to C when C ^ SFq. 
This demonstrates that the optimization of ZT can be 
realized by the enhancement of PF for a given Kpa- 

In Fig. 5, we consider square nanowires with a small 
cross section with Ly = Lz = 5 nm. To understand 
the dependence on QDA cross-sectional area, we show 
in Fig. 6 ZT as a function of temperature for various 
values of Lz with Ly fixed at 5 nm. We adopt F^ = 
F/i = F = 3ro, tc = 3.5ro, U( = 30ro and A = dOFo. 
In our model, we keep a large separation between QDAs 
along the z direction {Dg = 25 nm), this allow us to 
avoid the electron hopping effects and interdot Coulomb 
interactions between QDAs along that direction. Fig. 6 
shows the maximum ZT of QDA with Lz = 7 Ly can be 
larger than one. This is due to the enhancement of Ge 
with increased tunneling rates (F = SFo). The ZT value 
reduces significantly for the cases with large Lz (> 13Ly). 
This ZT reduction comes from the increased value of Kpa 
[see Fig. 1(b)]. If we reduce the separation, Dg between 
QDAs to increase the density of QDAs, the proximity 
effect arising from interdot Coulomb interactions between 
QDAs should be included. [33] In general, the proximity 
effect from interdot Coulomb interactions will suppress 
the Ge values. 

IV. SUMMARY AND CONCLUSIONS 

We have theoretically investigated the TE properties 
of QDAs embedded in nanowires connected to metallic 
electrodes in the Coulomb blockade regime. Both OCA 
and HEA calculations were examined. It was shown that 
HFA method gives unphysical results for the Seebeck co¬ 
efficient at low temperature when the QDs are nearly half 
hlled. However, it gives reasonable results at high tem¬ 
peratures when the QD energy levels are far above the 
Fermi level (low Filing condition). Under the condition 
Kpa S> Ke, the Optimization of ZT can be achieved by re¬ 
ducing Kpa and increasing the power factor PF = S^Ge, 
simultaneously. We found that for various design param¬ 
eters, ZT > 1 can be achieved. In our simulation, we 
assumed that the presence of QDs reduces the phonon 
conductivity in nanowires by a factor F), = 0.1. It is con¬ 
ceivable that a smart nanostructure design which reduces 
the value of Fg further can increase the maximum ZT 
even more. The maximum ZT value of a QDA nanowire 
with small cross section is found insensitive to tempera¬ 
ture variation (see Figs. 5(b) and 6). Such a feature is 
quite different from that of conventional thermoelectric 
materials. [1,2] The ability of QDAs to maintain ZT near 
its maximum value for a wide range of temperatures is 
an important feature in the realization of useful thermo¬ 
electric devices. [2] 
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FIG. 1: (a) Phonon thermal conductivity of a silicon nanowire 
as a function of temperature for the different cross sections. 
The square cross section is defined as A = Ly x Lz, where Ly 
and Lz are lateral sizes, and (b) phonon thermal conductivity 
of a GaAs nanowire with rectangular cross section A — Ly x 
Lz as a function of temperature for different Lz values at fixed 
Ly = 5 nm. 


FIG. 2: (a) Electrical conductance (Ge), (b) Seebeck coef¬ 
ficient (S), (c) electron thermal conductance (ace), and (d) 
power factor PF = S^Ge of triple QDs junction systems as a 
function of gate voltage for two different temperatures. Solid 
line and dot-dashed line are calculated by OCA and HFA, 
respectively. 
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FIG. 3: (a) Electrical conductance (Ge), (b) Seebeck coef¬ 
ficient (S), (c) electron thermal conductance (ke), and (d) 
power factor PF = S^Ge of triple QDs junction systems as a 
function of detuning energy A = Eo — Ef for different tem¬ 
peratures. The other physical parameters are the same as 
those of Fig. 2. 



FIG. 4: (a) Electrical conductance, (b) Seebeck coefficient 
and (c) {ZT)o of a QDA system with dot number N=10 in 
the absence of Uph as a function of temperature for different 
detuned energies. Other physical parameters are the same as 
those of Fig. 2. 
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FIG. 5: Figure of merit of a QDA system with dot number 
N=10 in the presence of Kph as a function of temperature. The 
curves of diagram (a) correspond to those of Fig. 4. Diagram 
(b) shows ZT for various electron hopping strengths. 



FIG. 6: {ZT) as a function of temperature for different Lz 
values. The physical parameters are tc = S.SFq and Tl = 
r_R = SFo. Solid line {Lz = 7Ly), dashed line {Lz = IZLy), 
dotted line {Lz = 79Ly) and dot-dashed line {Lz = 25Ly). 
Other physical parameters are the same as those of Fig. 5(b). 












